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ADMINISTRATIVE  INFORMATION 


The  Coast  Guard  Advanced  Marine  Vehicles  (AMY)  Project  is  under  the  direction 
of  the  Marine  Vehicle  Technology  Branch  (G-DMT-2)  of  the  Marine  Technology 
Division  (G-DMT)  in  the  Office  of  Research  and  Development  (G-D)  in 
Washington,  DC.  The  AMV  Project  (9207)  consists  of  three  major  elements; 
Operations  Research  (9207.1),  Ship  Test  and  Demonstration  (9207.2),  and 
Hydrodynamics  (9207.3).  The  Search  and  Rescue  Model  described  in  this  report 
is  part  of  the  Operations  Research  element.  It  was  developed  by  the  Marine 
Systems  Branch  of  the  Coast  Guard  Research  and  Development  Center  in  Groton, 
Connecticut.  The  Research  and  Development  Center  is  a  Headquarters  Unit 
reporting  to  the  Chief  of  the  Office  of  Research  and  Development  (G-D). 


INTRODUCTION 


The  U.S.  Coast  Guard  is  in  the  process  of  evaluating  Advanced  Marine 
Vehicles,  such  as  hydrofoils,  planing  boats.  Small  Waterplane  Area  Twin  Hulls 
(SWATH)  and  surface  effect  ships,  as  potential  Coast  Guard  cutters.  A 
necessary  part  of  this  process  is  to  quantify  the  operational  performance  of 
AMV's  and  conventional  displacement  craft  in  Coast  Guard  missions  so 
comparisons  can  be  made  between,  and  among,  these  different  types  of  vessels. 
Quantifying  performance  is  only  part  of  the  problem.  To  take  advantage  of 
their  unique  characteristics,  AMV's  may  be  operated  differently  than 
conventional  cutters.  It  is  also  necessary  to  understand  how  these  different 
methods  of  operation  may  affect  the  operational  performance. 

An  accepted  method  of  addressing  these  concerns  is  by  the  use  of  a  model. 
To  be  useful,  a  model  must  be  flexible  and  still  address  the  technical  issues 
of  the  problem.  In  the  case  of  Search  and  Rescue  (SAR),  the  model  must  be 
able  to  address  the  salient  aspects  of  a  typical  case  such  as  distance  to 
datum,  search,  survival  time,  and  weather,  but  more  importantly,  the  model 
must  produce  outputs  that  relate  directly  to  the  goal  of  the  SAR  program 
(saving  lives).  Furthermore,  the  outputs  should  be  quantifiable  and 
verifiable  in  practice.  The  goal  of  this  model  of  SAR  is  to  be  able  to 
quantify  the  lifesaving  performance  of  various  AMV  concepts  and  understand  how 
their  unique  characteristics  influence  mission  performance.  The  rationale 
for  quantifying  SAR  performance  for  a  vessel  is  presented  in  Reference  1. 

Several  requirements  were  placed  upon  the  development  of  the  SAR  Model. 
They  were: 

1.  be  supportable  with  available  or  measurable  data, 

2.  be  capable  of  distinguishing  between  craft  types,  and 

3.  be  able  to  accurately  reflect  important  mission  factors. 

The  philosophy  was  to  model  the  vessel  and  operation  at  as  high  (macro)  a 
level  as  possible  and  still  satisfy  the  requirements. 

At  this  point,  it  is  important  to  consider  a  fundamental  assumption  of  the 
model,  that  is,  life-threatening  SAR  cases  deteriorate  with  time.  If  the 
mariner  is  in  need  of  assistance,  the  situation  will  not  improve  on  its  own. 
The  boat  will  eventually  sink,  he  or  she  will  become  exhausted  and  hypothermia 
will  set  in.  At  some  time,  the  maximum  survival  time,  tx,  the  mariner  will 

no  longer  survive.  We  model  the  probability  that  the  mariner  is  alive  with  a 

monotonically  decreasing  function  of  time  which  becomes  zero  at  tx.  This 
curve  is  called  the  survival  function.  The  smaller  the  value  of  tx,  the 
more  critical  time  is  in  the  SAR  case.  It  follows  that  if  tx  is  very  large, 
then  time  is  not  critical  in  the  case. 

Although  this  concept  of  a  survival  function  is  logically  consistent,  it 
generates  a  problem  that  cannot  be  ignored.  That  is,  it  is  very  difficult  to 
accurately  define  the  curve.  It  is  also  difficult  to  distinguish  between  the 
capabilities  of  vessels  when  the  cases  are  not  sensitive  to  time.  Moderate  or 

low  severity  SAR  cases  typically  exhibit  this  type  of  behavior.  An  attempt  by 

Arrigan  (Reference  2)  to  determine  survival  curves  for  life-threatening  (high 
severity)  cases  from  the  Coast  Guard  SAR  data  base  was  thwarted  by  the  lack  of 
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reported  information  and  the  limited  number  of  cases  with  the  same  set  of 
controlled  conditions.  For  an  initial  approximation,  the  hypothermia  curve 
from  the  National  SAR  Manual  (CG-308)  can  be  used  to  define  the  survival 
curve.  This  is  only  applicable  for  people-in-the-water  cases.  Furthermore, 
temperature  is  restricted  to  only  one  of  two  values  (over  or  under  60°F). 
Relative  comparisons  of  vessels  will  be  sensitive  to  the  survival  function. 
Even  though  the  shape  of  the  curve  and  maximum  survival  time  are  not 
accurately  known,  we  can  still  use  the  concept  and  vary  the  value  of  tx  as  a 
model  parameter.  With  successive  runs  of  the  model,  the  user  can  determine  a 
relative  measure  of  SAR  effectiveness  for  the  vessels  of  interest. 

The  original  model  of  SAR  used  multiple  convolutions  to  arrive  at  a 
solution.  The  approach  presented  here  uses  a  combination  of  discrete  and 
continuous  distributions  to  represent  the  variables  in  a  SAR  case.  This 
required  the  use  of  Monte  Carlo  techniques  to  draw  samples  from  the 
distributions  since  multiple  convolution  of  the  mixed  probability  densities 
appeared  intractable.  However,  the  entire  model  is  of  such  a  size  that 
thousands  of  replications  can  be  run  at  minimal  cost.  The  problem  of  data 
availability  is  attacked  essentially  by  breaking  down  the  critical  variables 
into  more  elemental  (and  thus  more  easily  estimatable)  factors.  For  instance, 
search  time  is  the  result  of  datum  uncertainty,  craft  search  capability 
(search  width  and  search  speed),  and  drift  factors  affecting  the  original 
datum.  Mathematical  search  theory  is  then  used  to  recombine  these  variables 
to  yield  the  distribution  of  search  times. 


SAR  MODEL 


Assumptions 

Before  proceeding  to  a  description  of  the  SAR  model,  it  should  be 
emphasized  that  the  function  of  the  model  is  not  to  do  a  complete  analysis  of 
the  SAR  mission.  Its  purpose  is  only  to  compare,  in  a  relative  sense,  the 
lifesaving  capability  of  different  craft  under  the  same  set  of  conditions. 
This  comparison  is  based  on  the  distribution  of  the  relevant  factors  of  a  SAR 
case,  such  as  sea  state  and  the  distance  to  the  reported  datum,  and  vessel 
characteristics  such  as  speed  in  seaway,  sensor  sweep  width,  and  reliability. 
It  was  necessary  to  make  certain  assumptions  to  make  the  model  results  more 
useful.  Specifically: 

1.  Saving  a  life  is  the  primary  goal  in  SAR,  salvaging  property  is 
secondary.  If  a  vessel  is  a  good  lifesaver  it  is  assumed  to  be  a  good 
property  salvor. 

2.  Structural  details  of  specific  craft  types  are  not  modeled.  The 
lifesaving  function  can  be  performed  equally  well  from  any  vessel. 

3.  Design  will  account  for  all  appropriate  safety  aspects.  Good 
seamanship  is  assumed. 

Specific  design  requirements  can  be  built  into  any  of  the  AMV  concepts. 
For  example,  if  low  speed  maneuverability  is  a  requirement,  an  out-drive  can 
be  installed  on  any  type  of  vessel.  If  low  freeboard  or  a  specific  height  of 
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eye  are  required  in  SAR,  they  can  be  incorporated  in  the  design  of  the 
vessel.  The  fundamental  difference  in  vessel  types  for  purposes  of  the  model 
is  the  maximum  speed  of  the  vessel  and  the  speed  it  makes  in  a  seaway. 

With  this  background,  the  computational  flow  of  the  model  can  be  traced 

(Figure  la  and  lb).  For  consistency,  boxes  with  double  borders  are  used  for 

distributions  of  variables  which  are  drawn  at  random,  calculations  and  other 
deterministic  relationships  are  in  single-bordered  boxes,  and  data  passed  from 
one  model  section  to  another  are  in  circles. 

Starting  at  the  top  left  of  Figure  la,  a  sea  state,  SS,  is  drawn  from  the 

distribution  of  sea  states  and  passed  to  the  box  to  the  right  containing  the 

vessel  speed-sea  state  characteristic  curve  of  L.  The  speed  in  the  seaway,  V, 
is  then  calculated.  After  drawing  a  distance  to  datum,  D,  for  the  SAR  case, 
the  speed,  V,  and  distance,  D,  are  combined  to  produce  a  sample  transit  time 
to  datum,  t] .  Samples  are  then  drawn  from  the  three  remaining  distributions 
on  the  left  of  Figure  la,  i.e..  Time  to  get  Underway,  Incident  to  Notification 
Delay,  and  Incident  Time  of  Day.  The  24-hour  clock  time,  Tp,  is  then 

computed  as  well  as  the  elapsed  time,  to,  from  the  inception  of  the  case 
until  the  vessel  arrives  at  the  datum. 

The  probability  of  survival  versus  the  time  (from  incident)  curve  is  in  the 
upper  right  of  Figure  la.  This  curve  is  chosen,  as  are  all  the  probabilities 
in  the  model,  to  reflect  the  class  of  scenarios  being  considered.  The  time 
tx  is  an  estimate  of  the  maximum  survival  time,  given  the  environmental 
conditions  and  the  type  of  incidents  being  investigated  (e.g.,  persons  in  life 
raft,  in  the  water,  etc.).  The  search  time  is  computed  by  assuming  that  the 
search  will  be  terminated  at  the  next  sunset  after  the  maximum  survival  time 
has  elapsed  (given  that  the  search  has  been  unsuccessful  to  that  point). 

There  are  certainly  remarkable  cases  on  record  where  people  have  survived  for 
extended  periods  of  time  "against  all  odds,"  but  it  is  equally  true  that  every 
unsuccessful  search  is  stopped  at  some  point.  The  rule  used  here  to  stop  the 
search  was  chosen  as  a  reasonable  one,  and  is  consistent  with  the  overall  goal 
of  the  model,  that  is,  to  fairly  compare  various  craft  in  SAR. 

The  next  step  in  the  model  is  to  prepare  a  probability  of  detection  curve 
extending  up  to  the  time  the  search  is  terminated.  The  results  of  Koopman 
(Reference  3)  on  optimal  search  of  bivariate  (two  variables)  normal  target 
distributions  are  used,  along  with  a  splicing  technique  to  account  for 

sequential  changes  in  search  capability  in  daylight  and  darkness.  An  initial 
estimate  of  target  location  is  quantified  in  terms  of  the  model  input  CEP 

("circular  error  probable",  or  the  radius  of  the  smallest  circle  with  a  50% 
chance  of  including  the  actual  target  location).  This  circular  distribution 
is  modified  by  the  drift  effects  of  current  and  wind  (labeled  oc  and<rw  in 
the  lower  left  corner  of  Figure  lb)  to  produce  an  expanding  target 
distribution.  The  details  of  these  calculations  are  given  in  Appendix  A.  The 
craft's  search  capabilities  are  input  in  terms  of  search  speed  and  search 
width  for  the  two  different  regimes  (day/night).  All  of  this  information  is 
mathematically  combined  to  produce  probability  of  detection  (POD)  curves  for 
day  and  night.  These  curves,  in  conjunction  with  the  data  on  sunrise  and 
sunset  times,  the  time  of  search  initiation  and  maximum  search  time,  yield  a 
composite  cumulative  probability  of  detection  curve  ("cumulative"  box,  Figure 
lb)  for  the  case  being  evaluated.  The  time  at  which  the  search  will  be 
terminated  is  labeled  Tj. 


SEARCH  AND 
RESCUE  MODEL 


To  generate  a  sample 
the  interval  [0,  1.0] 
search  time,  then  it 
Otherwise,  the  inverse 
find. 


time  to  detection,  a  number  is  drawn  at  random  from 
.  If  this  number  is  larger  than  the  POD  at  maximum 
s  assumed  that  the  search  is  terminated  unsuccessfully. 
POD  function  is  evaluated,  producing  the  sample  time  to 


The  time  to  arrive  on  scene,  t$,  is  computed  by  adding  the  time  to 
datum,  tQ,  to  the  time  to  find  tp.  The  probability  of  survival  is 

determined  by  entering  the  survival  curve  and  reading  off  the  ordinate  value 
associated  with  t$. 

Finally,  a  random  draw  is  made  from  the  craft's  reliability  curve  to  see 
if  it  is  operational  for  the  time  required  to  effect  the  rescue.  This 
information  is  not  presently  used  to  affect  the  results  of  a  case.  No  delays 
for  repair  are  yet  included  in  the  model.  Only  a  1/0  indicator  (vessel  up  or 
down  for  t$  hours  of  operation)  is  logged.  The  various  times  and 

probabilities  for  this  replication  are  then  stored  in  memory  for  future 
post-processing,  and  the  entire  procedure  is  begun  again. 

Computer  Program 

A  FORTRAN  program  has  been  written  and  tested  to  implement  the  SAR  model. 

A  listing  is  included  in  Appendix  C.  It  is  presently  running  on  a  Digital 

Equipment  Corporation  VAX  11/780  located  at  the  Naval  Underwater  Systems 
Center  in  New  London,  Connecticut.  In  operation,  the  user  is  prompted  to 
select  from  a  variety  of  available  distributions  to  reflect  the  scenario  and 
operating  capabilities  under  investigation.  In  addition,  the  user  has  the 
option  of  generating  new  distributions  if  it  is  felt  that  those  provided  are 
unsuitable.  Such  inputs  as  wind,  current,  search  widths,  and  sunrise  and 
sunset  times  are  completely  user  determined.  A  text  file  containing  general 
information  and  operating  instructions  is  available  in  the  same  user  area. 
Since  the  computer  is  located  at  a  secure  naval  facility,  access  is  restricted 
to  authorized  personnel.  Individuals  wishing  to  run  the  program  should 
contact  the  Marine  Systems  Branch  at  the  Research  and  Development  Center  for 
specific  access  procedures. 

Post  Processing 

The  basic  procedure  in  the  data  reduction  of  the  outputs  of  Monte  Carlo 
models  consists  in  first  determining  a  quantity  to  be  estimated  (e.g.,  search 
time,  total  time  to  scene,  fraction  of  lives  saved),  averaging  it  over  some 
number  (n)  of  samples,  and  statistically  determining  a  confidence  interval 
containing  the  average  value.  Distributions  of  outputs  are  assumed  to  be 
asymptotically  normal.  Therefore,  the  confidence  interval  is  inversely 
proportional  to  the  square  root  of  the  sample  size.  If  the  confidence 
interval  is  too  large,  additional  samples  can  be  run.  Since  there  is  a 
relationship  O/vHi)  between  the  size  of  a  confidence  interval  and  the  number 
of  samples,  an  estimate  can  be  made  of  the  number  of  additional  samples 
necessary  to  yield  a  confidence  interval  of  the  desired  size.  The  model  can 
be  run  again  for  the  newly  computed  sample  size  to  produce  a  new  point 
estimate  and  confidence  interval. 


A  program  to  perform  the  statistical  processing  on  fraction  of  lives  saved 
is  currently  available  in  the  same  file  space  in  which  the  model  resides. 


Details  of  its  operation  and  user  information  for  the  SAR  model  can  be  found 
in  the  accompanying  text  file. 


Output  Format 

The  model  was  developed  to  compare  different  craft  in  various  operational 
scenarios.  Two  basic  styles  of  one-page  reports  were  designed  to  convey  the 
results  as  concisely  as  possible.  Figure  2  is  an  example  of  the  first  of 
these,  and  should  be  self-explanatory,  given  the  preceding  discussion.  The 
mission  performance  section  gives  the  bottom  line,  that  is  the  overall 

fraction  of  lives  saved.  The  remainder  of  cases  are  broken  out  into  the 
fraction  of  those  cases  where  the  search  resulted  in  finding  the  subject 
expired,  and  that  fraction  where  the  search  itself  was  unsuccessful  in 

locating  the  mariner.  The  overall  results  then  sum  to  1.0.  The  95% 
confidence  interval  is  that  obtained  by  the  usual  t-statistic,  i.e.,  the  model 
output  is  considered  normally  distributed  about  the  mean  value  of  fraction  of 
lives  saved. 

The  second  output  is  the  graphical  layout  of  Figure  3.  The  central  body 
of  the  graphic  depicts  the  scenario,  highlighting  the  craft  capabilities  on 
the  left  and  the  elliptical  search  area  on  the  right.  The  larger  rectangle 
and  the  inset  cross-hatched  rectangle  show  day  and  night  search  speeds 
(horizontal  dimension)  and  sensor  search  widths  (vertical  dimension).  The 
product  of  speed  and  sweep  width  is  the  correct  measure  of  search  capability, 

and  the  areas  of  the  rectangles  show  this.  The  average  distance  to  datum  and 

the  average  time  to  datum  are  near  the  bottom  of  the  graphic.  The  elliptical 
region  depicts  the  minimum  area  with  a  50%  chance  of  containing  the  target 
(elliptical  equivalent  of  CEP)  and  its  area  is  noted.  The  cross  is  annotated 
with  the  maximum  survival  time. 

Across  the  top  of  the  graphic  are: 

1.  Specific  information  on  the  craft,  including  speed  and  range 
capabilities. 

2.  A  clock  showing  times  of  light  and  darkness  (night  is  cross-hatched) 
along  with  two  radial  lines  representing  quartiles  of  the  time  of  day  of 
incident  distribution  (in  this  example,  between  approximately  1300  and  1700 
hours  50%  of  the  cases  occur).  The  time  of  the  peak  caseload  is  indicated  by 
the  arrow. 

3.  A  wide  panel  giving  indicators  of  the  environment:  water  temperature, 
sea  conditions  (quantified  by  the  percentage  of  waves  larger  than  three  feet), 
and  the  wind  and  current.  The  absolute  headings  of  wind  and  current  are  not 
important  since  the  size  and  shape  of  the  target  uncertainty  region  are 
affected  only  by  the  relative  angle  between  them.  (This  angle  is  depicted  in 
the  graphic.) 

4.  The  performance  of  the  craft  in  terms  of  fraction  of  lives  saved, 
along  with  a  similar  bar  graph  for  a  user  selected  baseline  craft.  The  bar  on 
the  left  represents  the  fraction  of  SAR  cases  with  a  successful  outcome  (i.e., 
life  saved),  the  bar  on  the  right  represents  that  fraction  where  the  search 
found  the  subject  after  death  occurred,  while  the  gap  between  the  two  bars 
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EXAMPLE  OF  RADC  SAR  MODEL  SUMMARY  OUTPUT* 


**  SAR  VESSEL 


TYPE:  BOEING  PW  VARIANT 


MAX  SPEED  IN  CALM  WATER 

48.0 

KTS 

MAX  SPEED  IN  SEA  STATE  3 

""  42.0 

KTS 

AVG  TIME  TO  GET  UNDERWAY 

0.10 

HRS 

AVG  TIME  TO  I NCAPAC I TAT I NG'TKTIURE 

650.0 

HRS 

DAYTIME  SEARCH  SWEEP  WIDTH 

30.0 

NMI 

NIGHTIME  SEARCH  SWEEP  WIDTH 

10.0 

NMI 

**  CONCEPT  OF  OPERATIONS 

SEARCH  UNTIL  NEXT  SUNSET  AFTER  MAXIMUM  SURVIVAL  TIME. 

OPTIMAL  SEARCH  OF  TWO-DIMENSIONAL  NORMAL  TARGET  DISTRIBUTION. 

SEARCH  SPEEDS: 

DAYTIME  38.0  KTS 

NIGHTIME  .  ~~  12.0  KTS 

**  SCENARIO 


ONE  SAR  CASE  AT  A  TIME;  NO  MULTI-UNIT  SEARCH  OR  PRIORITY  DIVERSIONS. 
*  CASELOAD: 


AVG  DELAY  BETWEEN  INCIDENT  AND  CG  NOTIFICATION 
TIME  OF  PEAK  CASELOAD 
ABSOLUTE  MAXIMUM  SURVlOT'TIHE 
*  OPERATING  AREA: 

AVG  DISTANCE  TO  DATUM 
DATUM  UNCERTAINTY  {50%“CETT 

CURRENT  “ . “  “ . ~ 

WIND 

RELATIVTANGLOETWEEN  CURRENT  AND  'OTJ - 

FRACTION  OF  WAVES  GREATER  THAN  3  FEET 

SUNRISE 

SUNSET 


0.13  HRS 
1500  HRS 
6.2  HRS 

65.0  NMI 
17.6  NMI 
3.5  KTS 
15.0  KTS 
110  DEG 
0.22 

0500  HRS 
1900  HRS 


**  RESULTS 


*  OPERATIONAL  PERFORMANCE 
AVG  TIME  TO  DATUM 

SEARCH  AREA  (50%  CEP)  UPON  ARRIVAL  AT  DATUM  _ 
AVG  SEARCH  TIME  _ 

*  MISSION  PERFORMANCE 

FRACTION  FOUND  BEFORE  SURVIVAL  LIMITS 
FRACTION  FOUND  AFTER  SURVIVAL  LIMITS 
FRACTION  NOT  FOUND  (PRESUMED  DEAD)  _ 

OVERALL _ _ 

*  STATISTICS 

95%  CONFIDENCE  INTERVAL  AROUND  FRACTION  ALIVE 
NUMBER  OF  MONTE  CARLO  SAMPLES  _ 


1.92  HRS 
1250  SQM I 
3.50  HRS 

ALIVE  DEAD 
0.73  0.04 

0.16 
0.07 

TT773 — 0727 


(0.70,  0.76) 
1000 


♦Representative;  not  actual  model  outputs 
FIGURE  2 


indicates  the  remaining  fraction  of  cases  where  the  search  was  unsuccessful 
(presumed  dead).  The  entire  scale  thus  has  length  1.0. 

Both  output  formats  try  to  compress  into  a  single  sheet  a  great  deal  of 
information.  The  actual  inputs  to  the  model  are  in  large  part  distributions. 
Thus  an  attempt  has  been  made  to  use  either  one  or  two  numerical  quantities, 
usually  average  values,  or  a  graphical  representation  to  summarize  these 
distributions. 


EXAI*>LE  RUNS 

The  outputs  of  the  SAR  model  depend  upon  the  distributions  that  make  up 
the  model.  When  representative  and  realistic  distributions  can  be  obtained, 
the  model  can  be  used  constructively.  The  model  will  be  useful  certainly  for 
it's  intended  purpose,  that  is  to  compare  the  effectiveness  in  SAR  of  various 
ship  designs.  It  may  also  be  useful  in  determining  the  effect  of  various  ship 
design  characteristics  upon  SAR  performance.  Although  the  model  was  not 
designed  to  explore  the  operational  aspects  of  SAR,  it  may  be  useful  in 
determining  how  parameters  such  as  distance,  notification  time,  and  darkness 
affect  operational  performance. 

An  example  of  how  the  model  might  be  used  to  see  the  effect  of  the  seaway 
upon  six  vessels'  SAR  performance  follows.  Remember,  although  there  are  not 
many  variables  in  the  model,  a  large  number  of  runs  may  be  required 

(combinatorial  explosion)  since  the  number  of  values  for  each  parameter  are 
multiplied  to  give  the  total  number  of  runs. 

The  distributions  shown  in  Figure  4  were  chosen.  The  incident  times  were 
spread  uniformly  over  the  24-hour  day,  sunrise  was  at  0600  and  sunset  at  1800 
hours.  Appendix  B  provides  an  explanation  of  the  technique  used  to 
parameterize  the  wave  height  distributions.  The  six  vessels  and  the 

associated  speed-wave  height  relationships  are  presented  in  Table  1.  Each 
model  run  consisted  of  2000  replications,  i.e,  individual  SAR  sorties.  The 
model  outputs  are  presented  in  Figure  5.  It  can  be  seen  that  as  the  seaway 
increases,  the  fraction  of  lives  saved  decreases  with  all  of  the  vessels.  The 
higher  values  for  the  hydrofoil  can  be  attributed  to  its  high  speed,  both  in 
calm  water  and  in  a  seaway.  This  follows  since  we  are  considering  only 
life-threatening  cases  as  reflected  by  the  survival  curve  (Figure  4). 

It  is  interesting  to  note  that  the  fraction  of  lives  saved  and  the 
speed-sea  state  relationship  correlate  very  positively  for  the  six  vessels. 

The  correlation  coefficients  are  presented  in  Table  2.  Many  people  have  had 

this  insight  into  SAR  performance  for  a  long  time.  What  makes  the  model 
useful  is  to  be  able  to  quantify  results  and  make  the  correlations  to  see  if 
there  is  a  relationship  between  variables.  This  can  help  one  to  understand 
more  fully  the  relationship  between  vessel  characteristics  and  mission 
performance.  Operational  choices,  such  as  speed  to  search,  and  sensors 
extending  sweepwidth,  may  also  be  analyzed  using  the  model. 


TABLE  1 


VESSELS  AND  SPEED-WAVE  HEIGHT  RELATIONSHIPS 


Wave 

Height 


Hydrofoil 


Planing 

Hull 


SWATH 


ft  34.4 

31.6 

27.6 

30.1 

33.3 

25.8 

27.6 

27.2 

32.3 

21.3 

27.3 

24.; 

31.4 

18.0 

26.8 

22.: 

30.6 

15.9 

25.9 

21.1 

29.8 

15.1 

24.8 

19.' 

5.0 

0.1* 

23.4 

17.' 

5.0 

0.1* 

21.7 

14.; 

Displacement  95-Footer 


FIGURE  5  -  Example  of  SAR  Program  Results 


The  model  results  can  be  used  to  compare  the  SAR  performance  of  various 
AMV  and  conventional  vessels.  An  agreed  upon  set  of  conditions,  or  range  of 
values,  would  be  entered  into  the  model  and  a  series  of  runs  would  be  made. 
The  ultimate  rank  ordering  of  different  vessels  requires  an  accepted  Measure 
of  Effectiveness  { MOE )  and  evaluation  procedure.  For  example,  the  MOE  might 
very  well  be  the  fraction  of  lives  saved,  but  this  might  be  combined  with  an 
MOE  from  non  life  threatening  cases  such  as  time  to  arrive  on  scene.  Any 
evaluation  requires  one  to  completely  and  carefully  spell  out  the  procedure. 
This  model  can  be  very  useful  in  that  process. 
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GLOSSARY 


Datum:  The  most  probable  location  of  the  search  object  corrected  for  drift  at 
any  particular  moment  during  the  mission. 

Distribution:  A  quantitative  description  of  the  relative  likelihoods  of 
various  possible  outcomes  of  an  experiment. 

Monte  Carlo  Simulation:  A  probabilistic  model  of  a  system  whose  solution  is 
obtained  by  repeated  sampling  from  appropriate  distributions,  rather 
than  by  closed  form  manipulation.  As  in  any  simulation,  results  are 
specified  by  statistical  estimates. 

Probability  of  Detection  (POD):  A  mathematical  function  which  reflects  the 
odds  of  success  or  failure  in  detecting  the  search  target,  under  the 
assumptions  of  a  given  search  pattern,  sensor  search  capabilities, 
etc.  POD  is  always  between  0  and  1,  e.g.,  a  POD  of  0.5  implies  a 
50-50  chance  of  detecting  the  target. 
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APPENDIX  A 


In  this  appendix,  the  method  used  to  predict  the  expansion  of  the  target 
uncertainty  region  is  presented.  Referring  to  Figure  A-l ,  point  Pi  is  the 
original  datum  and  Sg  is  the  original  CEP.  An  empirical  result  developed  by 
the  Oceanography  Branch,  USCG  R&D  Center,  gives  R,  the  new  CEP  around  new 
datum  ?2,  after  drifting  a  distance  D  as: 

R  =  1.1  V-3  D2  +  S2  +  S2 


where  Sy  is  the  self-localization  accuracy  of  the  vessel  trying  to  get  to 
datum. 

From  information  in  the  National  Search  and  Rescue  Manual,  it  was  determined 
that  the  important  contributions  to  the  drift  distance  D  are: 

1.  Leeway,  taken  to  be  in  the  direction  of  the  wind, 

2.  Wind  current  drift,  also  taken  to  be  in  the  direction  of  wind,  and 

3.  Sea  current  drift,  in  the  direction  of  current. 

Referring  to  Figure  A-2,  an  application  of  the  law  of  cosines  gives 

(DS)2  =  (CD)2  +  (L+WD)2  -  2 f CD ) ( L+WD)  cos  {  n  -  $  ) 

where  L  *  leeway  speed 

CD  =  sea  current  drift  speed 
WD  =  wind-driven  current  drift  speed 
0  =  angle  between  wind  and  current 
DS  =  drift  speed 

Then  D  =  DS  •  T  where  T  is  the  drift  time 

The  following  representative  numbers  were  taken  from  the  National  Search  and 
Rescue  Manual : 


L  =  (0.05)(wind  speed) 

WD  =  (.033) (wind  speed)  -  0.23 
or 

0.0 


(whichever  is  larger) 


CD  =  current  speed 

These  values  have  been  hard-coded  into  the  SAR  model,  but  are  easily  changed 
if  it  is  deemed  necessary. 

Example:  wind  speed  =  30  kts 
current  =  2  kts 
0=  600 


APPENDIX  B 


Historical  wave  height  data  are  often  presented  in  bar  chart  form,  as  in 
Reference  B-l.  These  charts  depict  the  relative  frequency  of  occurrence  of 
waves  in  a  number  of  wave  height  "bins,"  for  a  given  location  and  time  of 
year.  Comparison  of  these  bar  charts  with  plots  of  frequency  functions  of  the 
binomial  distribution  suggest  the  following  method  of  parameterizing  a 
collection  of  realistic  wave  height  distributions.  Select  an  appropriate 
number  of  bins  and  compute  the  values  of  the  binomial  frequency  function 

P{k)  *  (6  )•<  (1  -0  )n‘k  ,  k  =  0,1,. ...n-i, 


for  various  values  of  the  parameter  0(9  e  [0,1]).  Figure  B-l  shows  six 
such  distributions  each  with  8  bins  (n=7),  and  demonstrates  the  range  of  wave 
height  distributions  which  can  be  simulated.  In  fact  it  is  not  difficult  to 
find  actual  locations  whose  recorded  wave  height  distributions  are  very  close 
to  those  generated  by  this  method. 

The  SAR  model  was  run  with  a  sequence  of  such  distributions  with  the 
parameters  ranging  from  0.05  through  0.50  in  increments  of  0.05.  The  mean  of 
the  binomial  distribution  given  above  is  n  0,  and  completely  characterizes  the 
distribution.  Since  these  distributions  represent  recorded  significant 
(average  1/3  highest)  wave  heights,  the  term  "mean  significant  wave  height" 
was  used  to  specify  each  distribution. 

Reference: 

B-l. 
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Wind  and  Wave  Summaries  for  Selected  U.S.  Coast  Guard  Operating  Areas. 
Coast  Guard  RAO  Center  Report  No.  CG-D-11-83. 


BINOMIAL/WAVE  HEIGHT  DISTRIBUTION 

THETA  -  .10 

Parent  Occuranca 


t  9  5  7  •  11  19  19 

Wova  Haight  (Ft) 


BINOMIAL/WAVE  HEIGHT  DISTRIBUTION 

theta  -  jo 


1  3  9  7  9  11  13  is 

Wova  Haight  (ft) 


BINOMIAL/WAVE  HEIGHT  DISTRIBUTION 

THETA  -.30 


9  7  a  it 

Wova  Haight  (ft) 


BINOMIAL/WAVE  HEIGHT  DISTRIBUTION 

THETA  -  .15 


3  3  7  9  11 

Wova  Haight  (ft) 


BINOMIAL/WAVE  HEIGHT  DISTRIBUTION 

THETA  -.25 


3  7  9  11  13  19 

Wova  Haight  (ft) 


BINOMIAL/WAVE  HEIGHT  DISTRIBUTION 

THETA  -  .35 

Parent  Oecuranca 


FIGURE  B— 1 


3  7  I  11 

Wove  Haight  (ft) 


rvj  /nj  rr>  ro  rf  LT<  LT>  •£•  iT'  Is-  r-  H'  csi  f-l  ITi  f'-  <33 


$  ty  main. for 


PROGRAM  MAIN 


COMMON /ARRAY/ SEA i 8.2 » . 

I  PERCNT. 

3FEED<  8) . 

DUMMY 1. 

DTOD<  8 . 2  > . 
DATUM. 

UNDRWYi 10.2) . 
UNDER . 

XNOTFY < 9,2', 
DELAY . 

KING ID ( 24.2). 
FEAK, 

SURVIVf  0:13.2) 
EXPIRE. 

DWND3K ( 2 ) 

COMMON / GARBAGE / I SEED , 

FT  I  ME . 
CURRENT , 
WIND. 

I ANGLE. 

I PEAK. 

ISCRN. 

MMAX 


CHARACTER* 1  CHAR 
CHARACTER* 2 4  FILE1.FILE2 


COMMON / FACTS / S IGMAX . 

2  3IGMAY , 

3  WIDTH ( 2 ) . 

4  SSPEED!  2 ) 
COMMON /UN IT /I CRT, 

2  ITERM 

I CRT =6 
ITERM =6 


CALL  SARI 
CALL  SAR2 


'SEA  STATE  ARRAY 
'PERCENT  TIME  OVER  SS  5 
'.SPEED  VS  SEA  STATE 
! AVAILABLE 

*  DISTANCE  TO  DATUM  ARRAY 

'MEAN  DISTANCE  TO  DATUM 

'.TIME  TO  UNDERWAY  ARRAY  DIMENSION 

•MEAN  TIME  TO  UNDERWAY 

' INCIDENT  TO  NOTIFY  TIMES 

'MEAN  DELAY  TIME  BEFORE  NOTIFICATION 

'.  INCIDENT  TIME  OF  DAY  ARRAY 

'TIME  OF  FEAK  CASELOAD 

! SURVIVAL  PROB.  ARRAY 

'TIME  OF  CERTAIN  DEATH!'. 

! DAWN / DUSK  TIMES 
'RANDOM  NUMBER  SEED 
! MEAN  TIME  TO  FAILURE 
'.OCEAN  CURRENT  KTS . 

'.WIND  VELOCITY  KTS. 

'.ANGLE  BETWEEN  WIND  &  CURRENT 
'.INTEGER  VERSION  OF  PEAK 
) OUTPUT  TO  SCREEN  INDICATOR 
'.TOTAL  NUMBER  OF  MONTE  CARLO  REPS 


'SEARCH  SPACE  SIGMA 
'.SEARCH  SPACE  SIGMA 
! DAY/NIGHT  SEARCH  WIDTHS 
'DAY /NIGHT  SEARCH  SPEEDS 
(WRITE  UNIT  NUMBER 
'READ  UNIT  NUMBER 


'J  X»  4-  to  CD  j  X  X*  cn  4-  .A 


$  ty  sari. for 


SUBROUTINE  3AR1 


COMMON  ARRAY / SEA< 8.2) . 

1  P'ERCNT . 

2  SPEED (8). 

2  DUMMY 1 . 

3  DT0D( 8,2). 

3  DATUM, 

UNDRWY  < 10,2 ) . 
UNDER . 

XNOTFY (9,2), 
DELAY , 

XINCIDi 24,2) . 
PEAK, 

3URVIV( 0:13,2) . 
EXPIRE . 

DWNDSK ( 2) 

COMMON /GARBAGE/ ISEED . 

FTIME , 

CURRENT . 

WIND. 

I ANGLE . 

I PEAK. 

IS CRN , 

3  MMAX 

CHARACTER* 1  CHAR 
CHARACTER* 24  FILE1.FILE2 


COMMON  /  FACTS  /  SIGMAX . 

2  3IGMAY, 

2  CEP, 

?  WIDTH ( 2 ) , 

4  3 SPEED <2) 

COMMON /UNIT/ ICRT. 

2  ITERM 

I CRT =6 
I TERM =0 


USER  ENTERS  DISTRIBUTION  CHOIC 


260u  WRITE(  ICRT,*  )  CHOICES  FOR  SEA 
WRITE ( ICRT,* '  A.  MEAN  3IG  WAVE 
WRITE' ICRT.* )  B.  MEAN  SIG  WAVE 
WRITE' ICRT,* )  C.  MEAN  SIG  WAVE 


!  SEA  STATE  ARRAY 
(PERCENT  TIME  OVER  SS  3 
(SPEED  VS  SEA  STATE 
(AVAILABLE 

(DISTANCE  TO  DATUM  ARRAY 

1  MEAN  DISTANCE  TO  LATUM 

(TIME  TO  UNDERWAY  ARRAY  DIMENSION 

1  MEAN  TIME  TO  UNDERWAY 

'INCIDENT  TO  NOTIFY  TIMES 

(MEAN  DELAY  TIME  BEFORE  NOTIFICATION 

1  INCIDENT  TIME  OF  DAY  ARRAY 

'TIME  OF  PEAK  CASELOAD 

(SURVIVAL  PROB.  ARRAY 

(TIME  OF  CERTAIN  DEATH! ! 

'DAWN /DUSK  TIMES 
(RANDOM  NUMEER  SEED 
(MEAN  TIME  TO  FAILURE 
(OCEAN  CURRENT  KT5 . 

'WIND  VELOCITY  KTS . 

(ANGLE  BETWEEN  WIND  0  CURRENT 
(INTEGER  VERSION  OF  PEAK 
(OUTPUT  TO  SCREEN  INDICATOR 
'TOTAL  NUMBER  OF  MONTE  CARLO  REPS 


SEARCH  SPACE  SIGMA 
SEARCH  SPACE  SIGMA 
CIRC  ERROR  PROBABLE 
DAY 'NIGHT  SEARCH  WIDTHS 
DAY 'NIGHT  SEARCH  SPEEDS 
WRITE  UNIT  NUMBER 
READ  UNIT  NUMBER 


•TATE  ARE: 
HEIGHT- 2.0  FT 
HEIGHT- r.  '  FT 
HEIGHT =6.:  FT 


WRITE (I CRT,*) 'ENTER  A,  B.  OR  C 
READ( ITERM, 3601) CHAR 
FORMAT (Al) 

IF ( CHAR . NE.  A' . AND . CHAR . NE . 'B' . AND . CHAR . NE .  C' ) THEN 
WR ITE ( ITERM ,  * )  ' BAD  CHOICE' 

GO  TO  3600 
END  IF 

IF (CHAR. EQ.  A' (THEN 

OPEN ( UNIT* 2, NAME* ' STATE1 ' . FORM =  FORMATTED' ,TYPE=  OLD' ) 
END  IF 

IF ( CHAR . EQ . 'B’ )THEN 

OPEN(  UNIT=2  .NAME*  '  STATE2  '  . FORM=  '  FORMATTED  .TYPE*  OLD'  ) 
END  IF 

IF ( CHAR . EQ .  C' >THEN 

OPEN ( UN IT = 2, NAME =  ' STATE 3' . FORM =' FORMATTED ' .TYPE* ’OLD'  > 
END  IF 

READ< 2 . * ) SEA , PERCNT 
CLOSE (UNIT* 2) 

IPRCNT=INT( 100 . * PERCNT ) 


SPEED  VS  SEA  STATE  CHOICES 


WRITEi ICRT,* i ' SPEED  VS  SEA  STATE  iI.E.  SHIP  TYPE)  CHOICES  ARE 
WRITE( ICRT , * )  A.  HYDROFOIL 
WRITE) ICRT.* >'B.  PLANING  HULL 
WRITE( ICRT, * >  '  C.  SWATH' 

WRITEI ICRT ,  * ) ' D.  SES ' 

WRITE < ICRT,*  > ' E.  DISPLACEMENT' 

WRITE(ICRT,*>  E.  95  FOOT  WPB' 

WR ITE ( ICRT .  * )  '  ENTER  A,  B.  C,  D,  E,  OR  F' 

READ( ITERM, 360 1 ) CHAR 

IF ( CHAR . NE.  A' . AND . CHAR . NE .  B' . AND . CHAR . NE .  C 

.AND.  CHAR. NE. 'D' .AND. CHAR. NE.  E  .AND. CHAR. NE.  F' (THEN 
WR ITE ( ITERM ,  *  >  ' BAD  CHO ICE ' 

GO  TO  5800 
END  IF 

IF(CHAR.EQ.  A' )THEN 

0PEN(  UNIT*  2  .NAME*  '  SPEED  1  .FORM* ' FORMATTED' .TYPE*  OLD  > 

END  IF 

IF; CHAR. EQ.  B' )THEN 

OPEN) UNIT* 2. NAME*  SPEED 2  .FORM*  FORMATTED  .TYPE*  OLD  J 
END  IF 

IF' CHAR. EC.  C ' ) THEN 

:> PEN (  UNIT*:  .NAME*  SFEED3  .FORM*  FORMATTED  .TYFE*  OLD  ) 

END  IF 

IF ( CHAR . EQ •  D' 'THEN 

OPEN' UNIT* 2, NAME*  SFEED4  .FORM*  FORMATTED  .TYPE*  OLD  > 

END  IF 

IF (CHAR. EC.  E '  » THEN 

OFEN i UNIT* 2 , NAME  =  SPEEDS  .FORM*  FORMATTED  .TYPE*  OLD  • 

END  IF 


IF ( CHAR . EQ .  F  'THEN 

OPEN (UN IT =  2. NAME =  SPEEDS  .FORM=  FORMATTED ' .TYPE*  OLD' 
END  IF 

READ( 2  ' SPEED . DUMMY 1 

CLOSE i UNIT=2 ) 


DISTANCE  TO  DATUM  CHOICES 


WRITEf I CRT , * i ' 
WRITEl ICRT , *  1 
WRITE ( ICRT,* ) 
WRITEl ICRT.*) 
WRITE ( ICRT,* ) ' 
READ! ITERM, 360 
IF ( CHAR . NE.  A 
WRITE ( ITERM. 
GO  TO  ’’700 
END  IF 

IF  < CHAR . EQ .  A 
OPEN ( UNIT* 2 . 
END  IF 

IFlCHAR.EQ.  B' 
OPEN! UN IT* 2 , 
END  IF 

IF (CHAR. EQ.  C 
OPEN( UN IT =2 , 
END  IF 

READi  2  .  *  )  DTOD , 
CLOSE( UNIT=2 ) 


CHOICES  FOR  DISTANCE  TO  DATUM  DISTRIBUTION  ARE 

A.  PEAKED  AT  37.5  NM  OFFSHORE 

B.  PEAKED  AT  57.5  NM  OFFSHORE 

C.  PEAKED  AT  77.5  NM  OFFSHORE' 

ENTER  A,  B,  OR  C 

1 ) CHAR 

. AND. CHAR. NE. 'B' . AND . CHAR . NE .  C  >THEN 
*>  BAD  CHOICE' 


)THEN 

NAME= ’ DIST1  .FORM*' FORMATTED  .TYPE*  OLD'  > 
)THEN 

NAME= ■ DIST2 ' .FORM*  FORMATTED ' .TYPE*  OLD' ) 
i  THEN 

NAME=  DIST2 ' , FORM*  FORMATTED  .TYPE*  OLD  > 


DATUM 


TIME  TO  UNDERWAY  CHOICES 


WRITE',  ICRT.*-  CHOICES  FOP  TIME  TO  UNDERWAY  DISTRIBUTION  ARE: 
WRITE! ICRT,* )  A.  IMMEDIATE' 

WRITEl ICRT.* -  B.  20%  SPIKE  AT  0  THEN  TENT  CENTERED  AT  5  MIN. 
WRITE i ICRT,* )  C.  UNIFORM  ON  0  TO  12  MINUTES 
WRITEl ICRT.* »  ENTER  A.  6 .  OR  C 
READ( ITERM, 3601 ) CHAR 

IF ( CHAR . NE.  A ' . AND . CHAR . NE.  B  . AND . CHAP . NE.  C  ) THEN 
WRITEl  ITERM,  *  )  BAD'  CHOICE' 

GO  TO  0900 
END  IF 

IF ( CHAR .  EQ .  A  ) THEN 

0PEN( UNIT-2. NAME*  GETGOl  .FORM*  FORMATTED  .TYPE*  OLD  - 
END  IF 

IF ( CHAR . EQ . ' B' )THEN 

OPEN! UNIT* 2, NAME*  GETG02  .FORM*  FORMATTED  .TYPE*  OLD  ) 

END  IF 


noon  o  rj  o  o  o  o 


IF ( CHAR . EO .  C1 )THEN 

OPEN i UNIT- 2 . NAME- ' GETGO 3  . FORM-  FORMATTED  . TYPE-  OLD ' ) 
END  IF 

READ <  2  .  *  > UNDRWY , UNDER 
CLOSE (UN IT- 2) 


TIME  TO  NOTIFICATION  DISTRIBUTION 


2100  WRITE! ICRT.*>  CHOICES  FOR  TIME  TO  NOTIFICATION  DIST.  ARE: 
WRITE! ICRT,*)  A.  IMMEDIATE' 

WRITE! ICRT.+i'B.  UNIFORM  ON  0  TO  I  HOUR' 

WRITE! ICRT,*) 'C.  UNIFORM  ON  0-1  HR  THEN  RAMP  DOWN  TO  S  HRS' 
WRITE!  I.CRT.  +  !  ENTER  A.  B,  OR  C 
READ'!  ITERM, 3601) CHAR 

IF!  CHAR.  NE.  'A  .  AND .  CHAR  .  NE .  B'  .  AND .  CHAR. .  NE .  C'  >THEN 
WR ITE ( ITERM , * )  ' BAD  CHO I CE 
GO  TO  12100 
END  IF 

IF! CHAR. EQ.  A  ) THEN 

OPEN! UNIT- 2. NAME-  NOTIF1  .FORM-  FORMATTED' .TYPE-  OLD ' ) 

END  IF 

IF ( CHAR . EQ . ' B' 'THEN 

OPEN (UN IT- 2, NAME-  NOTIF2  .FORM-  FORMATTED' .TYPE-  OLD  j 
END  IF 

IF ( CHAR . EQ .  C ' ) THEN 

OPEN! UNIT-2. NAME-  N0TIF3' , FORM- ' FORMATTED ‘ .TYPE-  OLD' ) 

END  IF 

READ  <  2 , *  >  XNOTF Y . DELA  Y 
CLOSE! UNIT- 2) 


INCIDENT  TIME  OF  DAY  DISTRIBUTION 


14300  WRITE! ICRT.*) 'CHOICES  FOR  INCIDENT  TIME  OF  DAY  DIST.  ARE: 
WRITE! ICRT,*)  A.  UNIFORM  OVER  ENTIRE  DAY 
WRITE! ICRT, * >' B.  BROAD  PEAK  AT  1600  HRS' 

WRITE! ICRT,*) 'C.  BROAD  PEAK  AT  2000  HRS 
WRITE!  ICRT,  *  i  ENTER  A,  B.  OR  O 
READ! ITERM, 3601) CHAR 

IF!  CHAR.  NE.  A  .  AND. CHAR  .  NE.  B  '.  AND .  CHAR. .  NE.  C '  ) THEN 
WRITE! ITERM. *)  BAD  CHOICE 

GO  TO  14300 
END  IF 

IF' CHAR. EO-  A  'THEN 

OFENi UNIT-2, NAME-  INCTIME1 ' .FORM- ' FORMATTED' .TYPE-  OLD'  ) 
END  IF 

IF! CHAR. EQ.  B '  > THEN 

OPEN! UNIT-2 .NAME-  INCTIME2  .FORM-  FORMATTED' .TYPE-  OLD  ) 
END  IF 

IF' CHAR. EO-  C' )THEN 
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0PEN(UNIT=2 ,NAME= ' 
END  IF 

READ ( 2 , * ) XINC ID , PEAK 
CLOSE*  UNIT=2 ) 


'  INCTIME3  . F0RM  =  FORMATTED '  ,TYPE= ' OLD '  ) 


SURVIVAL  TIME  DISTRIBUTION 


C 

16600 


WRITE* I CRT , * ) 

WRITE* I CRT , *  > 

WRITE* ICRT,*)  B.  20  HOUR  MAX  SURVIVAL  TIME' 
WRITE* ICRT. *> 'C.  30  HOUR  MAX  SURVIVAL  TIME' 
WRITE* ICRT,*)  ENTER  A,  B,  OR  C 
READ* ITERM, 3601 ) CHAR 

IF*  CHAR. NE. 'A'  . AND . CHAR . NE .  B'  . AND. CHAR . NE. 
WRITE* ITERM,*)  BAD  CHOICE 
GO  TO  16600 
END  IF 

IF ( CHAR . EQ .  'A'  .'THEN 

OPEN ( UNIT= 2 , NAME= ' SURV1 ' , FORM= ' FORMATTED 
END  IF 

IF* CHAR. EQ. 'B' )THEN 

OPEN*  UNIT=2 ,NAME=  SURV2 ' ,FORM= ' FORMATTED' 
END  IF 

IF* CHAR. EQ. 'C' )THEN 

OPEN*  UNIT=2 ,NAME= ' SURV3 ' . FORM =' FORMATTED' 
END  IF 

READ*  2 , * ) SURVI V , EXPIRE 
CLOSE*  UNIT=  2 ) 


CHOICES  FOR  SURVIVAL  TIME  DISTRIBUTION  ARE: 

A.  10  HOUR  MAX  SURVIVAL  TIME' 

B.  20  HOUR  MAX  SURVIVAL  TIME' 

C.  30  HOUR  MAX  SURVIVAL  TIME' 

ENTER  A,  B.  OR  C 


(THEN 


SURVI' ,FORM=' FORMATTED' ,TYPE= 


SURV2 ' ,FORM=' FORMATTED' ,TYPE= 


SURV3 ' . FORM =  FORMATTED' ,TYPE= 


ALL  THE  OTHER  GARBAGE 


WRITE* ICRT,* )  ENTER  TIME  AT  DAWN  <24  HOUR  CLOCK) 

READ* ITERM , * , ERR  =  82) DWNDSK  *  1 > 

IDAWN=INT( DWNDSK* 1 ) ) 

IF* DWNDSK* 1 ) .LT. 0.0. OR. DWNDSK* 1 > .GT. 2400. 0) GO  TO  32 
DD=( DWNDSK* 1 •-< INT* DWNDSK* 1 ) /100. ; *100. ) i /60. 

DWNDSK* 1 )= INT* DWNDSK* 1 ) / 100. >  +DD 
GO  TO  85 

WRITE* ICRT,* ) '  HOURS  ENTERED  ARE  NOT  WITHIN  RANGE  OF  0000- 
GO  TO  80 

WRITE* ICRT,*)  ENTER  TIME  AT  DUSK  v  24  HOUR  CLOCK) 

READ* ITERM,*, ERR  =  86 (DWNDSK*  2) 

IDUSK= INT* DWNDSK*  2) ) 

IF*DWNDSK(2) .LT. 0.0. OR. DWNDSK* 2) .GT. 2400. 0)GO  TO  86 
DD= ( DWNDSK*  2 ) - ( INT*  DWNDSK*  2 ) / 100 . ) *  100 .  i )/60. 

DWNDSK* 2) =INT< DWNDSK* 2)/ 100. >+DD 
GO  TO  87 

WRITE* ICRT.* ) '  HOURS  ENTERED  ARE  NOT  WITHIN  RANGE  OF  0000 
GO  TO  85 


.GT. 2400.0) GO  TO 


67  CONTINUE 

£AAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAA*AAAAAAAAAAAAAAAAAAAAAAAAAAA 
CEP =.l* DATUM 

C^^^^^^^^^^^^^^^^^^AAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAaaAAAAAAAAAAAAAAAAAAAAAAAAA 

90  WRITE( ICRT,*) 'ENTER  DAYTIME  SEARCH  WIDTH  (MILES)' 

READ ( ITERM , * , ERR  =  92)  WIDTH(l) 

GO  TO  94 

92  WRITE( ICRT.*) 'UNSATISFACTORY  RESPONSE;  TRY  AGAIN' 

GO  TO  90 

94  CONTINUE 

95  WRITE ( ICRT.*) 'ENTER  NIGHT  SEARCH  WIDTH  (MILES)' 

READ< ITERM, *.ERR=97 )  WIDTH<2) 

GO  TO  99 

97  WRITE< ICRT.*) 'UNSATISFACTORY  RESPONSE;  TRY  AGAIN ' 

GO  TO  95 

99  CONTINUE 

AAAAAAAAAAAAAAAAAAAAAAAAAAAAA^AAAAAAA^AAAAAAAAA'rfrAAAA^AAAAAAyyAAAAAAAAAAAAAA^A 

100  WRITEf ICRT,*) -ENTER  A  RANDOM  NUMBER  GENERATOR  SEED  (INTEGER  1-1000) 
READ( ITERM, *.ERR= 102)  ISEED 

GO  TO  104 

102  WRITE< ICRT.*)  UNSATISFACTORY  RESPONSE;  TRY  AGAIN 

GO  TO  100 

104  CONTINUE 

CAAAAAAAAAAAAAAAAAAAAAaAAAAAAAAAAAAAAaAAAAAAAAAAAAAAAAAAAAAAAAAA*AAAAAAAAAAA 

105  WRITEf ICRT.*) 'ENTER  MEAN  TIME  BETWEEN  FAILURE' 

READ < ITERM . * , ERR  =1051)  FTIME 

GO  TO  1052 

1051  WRITE< ICRT.*) ' UNSAT;  TRY  AGAIN' 

GO  TO  105 

1052  CONTINUE 

(^VAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAaAAAAAAAAAAAAAaAAAAAAAAAAAAAAAAAAAAAAAAAA 

c 

107  WRITE( ICRT,*)  ENTER  CURRENT  ( KTS ) ' 

READ( ITERM, *,ERR= 107 1JCURENT 
GO  TO  1072 

1071  WRITE( ICRT,*) 'UNSAT' 

GO  TO  107 

1072  CONTINUE 


108  WRITE ( ICRT, * ) ' ENTER  WIND  SPEED  i KTS ) 

READ( ITERM, *,ERR= 1081) WIND 
GO  TO  1082 

1081  WRITE( ICRT,*) ' UNSAT' 

GO  TO  108 

1082  CONTINUE 


109  WRITEt ICRT, *  > ' ENTER  ANGLE  BETWEEN  WIND  AND  CURRENT  (DEG) 
READ ( ITERM ,*, ERR = 1 0 9 1 ) I ANGLE 
GO  TO  1092 

1091  WRITE ( ICRT , * ) ' UNSAT ' 

GO  TO  109 
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1092  CONTINUE 

C 

c 

IPEAK=INT( PEAK ) 

C106  WRITE* ICRT,* ) ' ENTER  NUMBER  OF  SAMPLES  PER  REPLICATION' 

C  READ*  ITERM,*  ,£RR=1061 )  ISMPLS 

C  GO  TO  1062 

Cl 061  WRITE* ICRT.*)  UNSAT;  TRY  AGAIN' 

C  GO  TO  106 

Cl 06 2  CONTINUE 

110  WRITE* I CRT.*) 'ENTER  NUMBER  OF  REPLICATIONS' 

READ* ITERM. * ,ERR=1101 > IREPS 

GO  TO  1102 

1101  WRITE* ICRT,* > 'UNSAT;  TRY  AGAIN' 

GO  TO  110 

1102  CONTINUE 
MMAX= IREPS 

(^AAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAA 

WRITE* ICRT,*) ' ENTER  NAME  OF  FILE  TO  HOLD  USER  CHOICES' 

READ* ITERM, 111 >FILE1 

111  FORMAT*  A24 ) 

0PEN(NAME=FILE1 ,TYPE= ' NEW' ,FORM= ' FORMATTED' ,UNIT=3) 

WRITE* ICRT,*)  ENTER  NAME  OF  FILE  FOR  PROGRAM  OUTPUT' 

READ ( ITERM , 111 ) FILE2 

C  OPEN* NAME=FILE2.TYPE=' NEW' ,FORM= ' FORMATTED' .UNIT=2) 

OPEN*  UNIT=4 ,NAME= ' I ALIVE '  , FORM =' FORMATTED' ,TYPE='NEW'  ) 

C  WRITE* 2,*) ISMPLS, IREPS 

WRITE* 3,*) SPEED* 1) ,SPEED<3> .UNDER 
WR ITE < 3 , * ) FTIME , WIDTH ( 1 ) , WIDTH ! 2 ) 

WRITE ( 3 ,*) DELAY , IPEAK . EXP IRE 

WR ITE ( 3 ,*) DATUM , CEP . CURENT , WIND , I ANGLE 

WRITE ( 3 , * ) IPRCNT , IDAWN , IDUSK 

WRITE*  3 , * ) ISMPLS , IREPS 

WRITE*  4 . * )MMAX 

115  WRITE* ICRT,*) 'DO  YOU  WANT  OUTPUT  ON  SCREEN?  ( Y/ N ) 

READ* ITERM, 3601) CHAR 

IF ( CHAR . NE.  Y' . AND . CHAR . NE .  N'  > THEN 

WRITE* ITERM,* )  HUH?' 

GO  TO  115 
END  IF 

IF* CHAR . EQ .  Y'  )THEN 

ISCRN=1 

END  IF 

IF* CHAR. EQ.  N' )THEN 

ISCRN=0 

END  IF 

IF ( ISCRN  .EQ.  1)  THEN 
WRITE* ICRT, 400) 

END  IF 

400  FORMAT*  '1') 

RETURN 

END 
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~-i  ri  oj  rn  -d*  IP  IP 


$  ty  sar2.for 


Copy  av  ’;1  j 

permit  felly  RyifJie 


D"'TC  rig0  5  no* 
t‘-p..,d  action 


SUBROUTINE  SAR2 


b 

7 

7 

8 


3 

4 

5 

6 

7 

8 


COMMON /ARRAY /SEA (8, 2) . 

PERCNT , 

SPEED ( 8 ) . 

DUMMY 1 , 

DT0D( 8,2) , 
DATUM, 

UNDRWY ( 10.2) , 
UNDER , 

XNOTFY (9,2), 
DELAY , 

XINCIDt  24,2) . 
PEAK 

SURVIV (0:13,2), 
EXPIRE, 

DWNDSK ( 2 ) 

COMMON/ GARBAGE/ I SEED , 

F"IME, 

CL  RRENT . 

WIND. 

I ANGLE . 

IPEAK. 

I 3 CRN. 

MMAX 


!  SEA  STATE  ARRAY 
'.PERCENT  TIME  OVER  5  FEET 
!  SPEED  VS  SEA  STATE 
(AVAILABLE 

(DISTANCE  TO  DATUM  ARRAY 

(MEAN  DISTANCE  TO  DATUM 

(TIME  TO  UNDERWAY  ARRAY  DIMENSION 

(MEAN  TIME  TO  UNDERWAY 

(INCIDENT  TO  NOTIFY  TIMES 

(MEAN  DELAY  TIME  BEFORE  NOTIFICATION 

! INCIDENT  TIME  OF  DAY  ARRAY 

(TIME  OF  PEAK  CASELOAD 

(SURVIVAL  PROB.  ARRAY 

(TIME  OF  CERTAIN  DEATH!! 

(DAWN /DUSK  TIMES 
(RANDOM  NUMBER  SEED 
(MEAN  TIME  TO  FAILURE 
(OCEAN  CURRENT  KTS . 

(WIND  VELOCITY  KTS. 

(ANGLE  BETWEEN  WIND  6.  CURRENT 
(INTEGER  VERSION  OF  PEAK 
!  OUTPUT  TO  SCREEN  INDICATOR 
’.TOTAL  NUMBER  OF  MONTE  CARLO  REPS 


r 


CHARACTER*!  CHAR 
CHARACTER  *  2  4  FILE1.FILE2 


2 

3 

4 

5 
b 

7 

6 
? 

3 

7 

b 

b 

b 

s 

4 
4 


COMMON/ INFO /SSTATE, 
V, 

DIST, 
TTRANS , 
TUNDER , 
TNOTFY , 
TCINC , 
TDATUM , 
TCDATM , 
TEXPIR, 
I  UP. 
ISERCH , 
TFIND , 

I CESS. 
TSCENE , 
TCCENE , 
FALIVE. 


'SAMFLE  SEA  STATE 
(CORRESP.  SAMPLE  SPEED 
(SAMPLE  DISTANCE  TO  DATUM 
(CORRESP.  TRANSIT  TIME 
’SAMPLE  UNDERWAY  DELAY  TIME 
(SAMPLE  NOTIFICATION  TIME 
’SAMPLE  CLOCK  TIME  OF  INCIDENT 
(CORRESP.  TIME  TO  DATUM 
(CLOCK  TIME  TO  DATUM 
(TIME  TO  DEATH  W/  PROB  I 
■SHIP  UP  ON  SCENE:  1=YES 
(SEARCH  REO'D:  1=YES 
(SAMPLE  TIME  TO  FIND 
(SUCCESSFUL  SEARCH:  1=YES 
’RESULTANT  TIME  TO  SCENE 
’.CLOCK  TIME  ON  SCENE 
(RESULTANT  PROB  OF  SURVIVAL 


fjOOOuO 


3  IALIVE, 

2  IQUIT 

COMMON / FACTS / S IGMAX , 

2  3IGMAY. 

2  CEP , 

3  WIDTH ( 2 ) , 

4  SSPEED( 2 ) 
COMMON /UN IT' ICRT , 

2  ITERM 


'.SURVIVAL  INDICATOR:  1=  ALIVE 

! 1 =3EARCH  BROKEN  OFF  UNSUCCESSFULLY 

! SEARCH  SPACE  SIGMA 

'.SEARCH  SPACE  SIGMA 

! CIRC  ERROR  PROBABLE 

'.DAY /NIGHT  SEARCH  WIDTHS 

'.DAY /NIGHT  SEARCH  SPEEDS 

'.WRITE  UNIT  NUMBER 

'.READ  UNIT  NUMBER 


TEXPIR  =  SURVIV(  13.1) 

ANGLE=FLOAT( I ANGLE) 

EEWAY=. 05+WIND 
WDRIFT=.033*WIND-. 23 
IF ( WDRIFT. LT. 0 . )WDRIFT=0. 

DS  Q  =  (  EEWAY+WDR I  FT  >  *  *  2  +C  UP.ENT *  *  2 

DSQ=DSQ-2*CURENT*(EEWAY+WDRIFT>*COSi  3.  14 15:- -ANGLE) 

IF ( DSQ . L£. 0 )THEN 

DSQ  =  0 

GO  TO  26 

END  IF 

DRIFT=SQRT( DSQ ) 

26  CONTINUE 

DRIFT= .  5*DRIFI>TEXPIR 
R  =  1 . l+SQRT) . 3*DRIFT**2+CEP**  2  +  9 ) 

3IGMAX= . 849322*R 
SIGMAY=SIGMAX 


START  THE  MAIN  REPLICATION  LOOF 
LOOP  INDEX  IS  MAIN' 


MAIN= 1 

999  CONTINUE 

C  WRITE) 6,*  >  MAIN=  .MAIN 

IQUIT=0 

3  STATE  =  S  AMPLE ( SEA . 8 . 1 SEED > 

C  WRITE) 6.* i  WAVE  HEIGHT =  .SSTATE 

V = SPEED i S STATE) 

C  WRITE) 6 , * )  V=' .V 

SSPEED) 1 ) =V 
SSPEED)  2 )  =V 

D I  ST  =  SAMPLE  (  DTOD  .8,1  SEEL' ) 

C  WRITE) 6,*) 'DIST=' ,DIST 

TTRANS=DIST/ V 

C  WRITE) 6 , * ) '  TTRANS=  .TTRANS 

TUNDER *S AMPLE ( UNDRWY. 10 . I SEED ) 


n  o o o  n  non o non 


WRITE ( 6 ,  *  ) ' TUNDER=  .TUNDER 
TNOTFY = SAMPLE ( XNOTFY , 9 . I SEED ) 
WR ITE ( 6 , *  > ' TNOTFY  =  .TNOTFY 
TCINC=SAMPL£(XINCID.24.ISEED) 
WRITE( 6 ,  *  >  TCINC= ' .TCINC 
TDATUM=TTRANS+TUNDER+TNOTFY 
TCDATM =TC  INC  -t-TDATUM 


C  WRITE <  6  ,  *  >  ' TEXPIR= '  . TEXPIR 

TMAX = TEXP I R - TDATUM 
WRITE( 6 , * ) ' TMAX=  .TMAX 


DETERMINE  TOTAL  SEARCH  TIME 


TNEW=AMOD( TCDATM. 24. ) 

IF ( TNEW  .LT.  DWNDSK( 2 )  -AND. 

2  TNEW  ,GT.  DWNDSK (  1 )  )  THEN  '.LIGHT  AT  DATUM 

IF ( DWNDSK ( 2  > -TNEW  . GT.  TMAX)  THEN  IDEAD  BEFORE  DUSK 

TSERCH= DWNDSK ( 2 ) -TNEW 
ELSE 

TSERCH = DWNDSK < 2 ) -TNEW+24 .  ! SEARCH  EXTRA  DAY 

END  IF 

ELSE  IDARK  AT  DATUM 

IF  (  (  24  .  -  i  TNEW-DWNDSK  (  2  )  >  )  .GT.  TMAX)  THEN  '.DEAD  BFOR  NEXT  DSK 
TSERCH=24. - i TNEW -DWNDSK! 2) ) 

ELSE  ! DEAD  AFTER  NEXT  DUSK 

TSERCH=48 . - (TNEW-DWNDSK i 2  > ) 

END  IF 
END  IF 

WRITE(b,*> ' TSERCH= ' .TSERCH 


ICESS=1 

CALL  SEARCH ( TNEW . TSERCH . DWNDSK , I SEED , TF I ND ) 
C  WRITE(6.*) 'BACK  IN  SAR ,TFIND= ' .TFIND 

IF( TFIND  .EQ.  -1.0) THEN 
C  WRITE( 6 , * ) ' DID  NOT  FIND' 

ISERCH=1 
TFIND=TSERCH 
ICESS=0 
TSCENE=999 . 

TCCENE=999 . 

UPTEST=TDATUM+TSERCH 
PALIVE=0 . 

IALIVE=0 

IQUIT=1 

EVEN = RAN (I SEED) 

GO  TO  2000 


C-ll 


END  IF 

IF(TFIND  .LT.  .25)  THEN 
ISERCH=0 
ELSE 

ISERCH=1 
END  IF 

TSCENE =TDATUM+TF IND 
TCCENE=TCDATM+TFIND 
UPTEST=TSC£NE 
J=1 

IF (TSCENE  .LT.  SURVIV(J,1>>  THEN 
PALIVE=SURVIV( J- 1 , 2 ) 

ELSE 

J=J+1 

IF( J  .GT.  13)  THEN 

WRITE( ICRT,*> 'TIME  TO  FIND  WAS  TOO  LONG;TFIND=  .TFIND 
PALIVE=0 . 

GO  TO  1001 
END  IF 
GO  TO  1000 
CONTINUE 
END  IF 

WRITE<6,*)  PALIVE= ' ,PALIVE 
TEST1=RAN( ISEED) 

IF(TEST1  .LT.  P ALIVE)  THEN 
IALIVE=1 
ELSE 

IALIVE=0 
END  IF 
CONTINUE 
RLAM-l./FTIME 
REL=EXP( -RLAM*UPTEST) 

TEST = RAN ( ISEED) 

IF (TEST  .GT,  REL)  THEN 
IUP=0 
ELSE 
IUP=1 
END  IF 

IF ( ISCRN  ,E0.  1)  THEN 
CALL  SCREEN (MAIN) 

END  IF 
CALL  FILE 
MAIN=MAIN+1 

IF (MAIN  .LE.  MMAX )  THEN 
GO  TO  999 
ELSE 

WRITE ( ICRT , *  )  'FINISHED' 

END  IF 

CL0SE( UNIT=2 ) 

RETURN 
END 


onn  o  n  o  onn o  o o 


SUBROUTINE  TO  WRITE  SAR  MOE  OUTPUTS  TO  THE  SCREEN 
MAIN  IS  THE  REPLICATION  LOOP  COUNTER  IN  PROGRAM  RWSAR 
THE  HEADINGS  ARE  REPEATED  EVERY  20  LINES 


SUBROUTINE  SCREEN)  MAIN ) 

COMMON / INFO / S STATE , V , DIST , TTRANS , TUNDER , 

2  TNOTFY , TC INC , TDATUM , TCDATM , 

3  TEXFIR , IUP . ISERCH . TFIND , ICESS , TSCENE , 

4  TCCENE.PALIVE , I ALIVE 
COMMON / UNIT / ICRT . ITERM 


IF(MOD( MAIN, 20 )  .EQ.  1)  THEN 

WRITE ( ICRT,*  )  '  I - CASE  DATA - I  | - OPERATIO 

2NS - I  | ---PERFORMANCE--- I ’ 

WRITE) ICRT,*)  CLOCK  TIME  DISTANCE  TIME  TO  TIME  TO  TRANSIT 
2TIME  SPENT  SUCCESSFUL  PROB  OF' 

WRITE)  ICRT,*) 'OF  INCIDENT  TO  DATUM  NOTIFY  UNDERWAY  TIME 

SEARCHING  SEARCH  SURVIVAL' 

END  IF 


WR ITE ( ICRT , 1 0  0 ) TC INC , DI ST , TNOTFY , TUNDER .  TTRANS , 

TFIND. ICESS. PALIVE 

FORMAT ( 3X , F6 . 2 , 4X , F6 . 2 , 3X , F6 . 2 , 2X , F6 . 2 , 3X , F6 . 2 , 4X , F6 . 2 , 3X , I 6 , 
4X.F6. 2) 

RETURN 

END 


2 


2 

100 


SUBROUTINE  SEARCH (TD, 

2  TS, 

3  TIME, 

4  I SEED, 

5  TTF ) 
DIMENSION  TPOD) 0:5)  , 

3  POD) 0:5). 

4  TT( 2 ) , 

5  LIGHT) 5). 

6  DTIME) 6,2), 

7  TIME) 2) 

C  WRITE (6,*)  NOW  IN  SUBROUTINE 


'.TIME  OF  DAY  AT  DATUM 
(TOTAL  ALLOWABLE  SEARCH  TIME 
(DAWN/DUSK  TIMES  (24  HR  CLOCK ) 
(RANDOM  NUMBER  SEED 
(TIME  TO  FIND  s  RETURNED ) 

(TIMES  OF  SUCC.  PHASES  OF  SEARCH 
(POD  AFTER  SUCC.  PHASES 
(LENGTH  OF  DAY /NIGHT 
! DAY  =  1 ,NIGHT=2  INDICATOR 
(START/STOP  TIMES  FOR  DENSE)*). 
(DAWN/ DUSK  TIMES  (24  HR  CLOCK) 
SEARCH' 
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WRITE ( 6 , * )  TD= ' ,  TD , '  TS  =  ,TS 


IF  TIME  AT  DATUM  IS  DURING  DAYLIGHT.  ISTART=1 


IF(TD  .LT.  TIME( 2 )  .AND.  TD  . 

I START =1 
ELSE 

I START =2 
END  IF 

WRITE (6,*) ' I START = ' .ISTART 


TPOD( 0 ) =0 . 

I=MOD( ISTART. 2 )+l 
J  =  0 

T=TIME( I) -TD 
IF( T  .LT.  0. )  THEN 
T=T+-24 . 

END  IF 

TT( 1 ) =TIME( 2 ) -TIME( 1 ) 

TT( 2 ) =24 . -TIME( 2 ) +TIME( 1 ) 

WR ITE  ( 6 ,  * )  '  LENGTH  OF  DAY/NIGI 


.  TIME( 1) 'THEN 


'.FOR  LOOPING  CONVIENIENCE 
(INITIALISE  DAY/NIGHT  COUNTER 
! INITIALIZE  TPOD  SUBSCRIPT 
(TIME  REMAINING  IN  DAY/NIGHT 


(LENGTH  OF  DAY 
(LENGTH  OF  NIGHT 
=' .TT( 1 ) ,TT( 2 ) 


CONTINUE 

IF(T  .GE.  TS)  GO  TO  22 


J  =  J+1 
TPOD( J)=T 
T=T+TT( I) 

I=MOD( 1 ,2  >  +  l 

LIGHT(J)=I 

WRITE ( 6 , * ) ' J=' ,J 

WRITE* 6,*) 'TPOD'  ,J,  ' =  '  ,TPOD( J) 

WRITE ( 6 , * ) 'NEXT  T= ' ,T 


(INCREMENT  TPOD  SUBSCRIPT 
(SET  TPOD 

(INCR.  T  BY  NEXT  DAY / NIGHT  LENGTH 
(TOGGLE  DAY /NIGHT  COUNTER 
(SET  DAY/NIGHT  INDICATOR 


WRITE* 6,* )  NEW  DAY/NIGHT  INDICATOR= ' , I 

GO  TO  21 

CONTINUE 


J=J  +  1 

I=MOD( I , 2 ) +1 
LIGHT* J ) =1 

TPOD( J ) =TS  'FINAL  TIME  IS  TOT.  ALLOW.  TIME 

WRITE* 6.* )  FINAL  TIME  =  .TS,' FINAL  DAY /NIOHT=  , I 


ooooo  non  *—  »-■  nno  noon  non 


C 


P0D( 0 ) =0 

DTIME( 1,1) =TP0D( 0 ) 
DTIME! 1,2) =TP0D< 1 ) 


!  FIRST  DENSE_TIMES  ARE  EASY 


DO  100  I«1.J 

IF(  I  . EQ.  1)  GO  TO  150 

DTIME( 1 , 1 )  = 

2  CINVl LIGHT ( I ) ,P0D< I- 1 ) , DTIME! 1-1.2) ) 

DTIME( I , 2 ) = 

2  DTIME (1,1) +TP0D ( I ) -TPOD ( I - 1 ) 

150  CONTINUE 

POD! I) =P0D< I-l)+ 

2  CUMULI  LIGHT (  I )  , DTIME ( 1,2)  )  - 

3  CUMULI  LIGHT! I) .DTIME! 1,1) ) 

WRITE! 6 , * ) '1=  .1 

WRITE! 6 , * ) 'DTIME! 1,1)=' , DTIME! 1,1) 

WRITE! 6,*) 'DTIME! 1,2)=' , DTIME! 1,2) 

WRITE! 6.*) ' POD! I )  =  .POD! I) 

100  CONTINUE 


U=RAN< ISEED ) 

IF  (U  .GT.  POD! J ) )  THEN 
TTF = - 1 . 

RETURN 
ELSE 
1  =  1 

02  CONTINUE 

IF! POD! I )  .GT.  U)  GO  TO  103 
1  =  1  +  1 
GO  TO  102 
03  CONTINUE 

TTF=TP0Dl 1-1 ) +CINV! LIGHT! I ) ,U, DTIME! 1,2) > -DTIME! 1,1) 
END  IF 


RETURN 

END 


FUNCTION  TO  COMPUTE  THE  CUMULATIVE  POD  AT  A  GIVEN  TIME. 

THE  POD  FUNCTION  IS  THAT  GIVEN  BY  K00PMAN  FOR  EXPONENTIAL 
SENSOR  CAPABILITY  AND  OPTIMAL  ALLOCATION  OF  SEARCH  EFFORT. 


C-15 


FUNCTION  CUMUL( I .TIME)  ! I =1  MEANS  DAY ; = 2 MEANS  NIGHT 

COMMON/ FACTS /SIGMAX,  'SEARCH  SPACE  SIGMA  iX) 

SIGMAY,  '.SEARCH  SPACE  SIGMA  (Y) 

CEP,  ! CIRC  ERROR  PROBABLE 

WIDTH! 2),  ! DAY/NIGHT  SEARCH  WIDTHS 

SPEED!  2)  '.DAY /NIGHT  SEARCH  SPEEDS 


PI=3. 14159 

C=WIDTH( I)*SPEED( I) / ( SIGMAX*SIGMAY*PI ) 
ARG=SQRT( C*TIME) 

CUMUL=1. -( l.+ARG)*(EXP< -ARG) ) 

RETURN 

END 


FUNCTION  TO  COMPUTE  THE  INVERSE  OF  THE  SPLICED  CUMULATIVE 
POD  OF  THE  SEARCH.  AS  THE  POD  FUNCTION  IS  TRANCENDENTAL . 
AN  ITERATIVE  TECHNIQUE  IS  USED. 

I =DAY/ NIGHT  INDICATOR ; 1 =DAY , 2=NIGHT 
POD=CUMULATIVE  POD  ON  DISTRIBUTION  CUMULI  I) 

TSTART= STARTING  TIME  FOR  ITERATION 


FUNCTION  CINV< I , POD , TSTART ) 


COMMON / FACTS/S IGMAX,SIGMAY, CEP, WIDTH! 2) , SPEED! 2) 
F( T) =1 . -POD- ( 1 . +SQRT! C*T) ) *EXP< -SORT! C*T) ) 


PI*3. 14159 

C= WIDTH!  I )  ’‘'SPEED!  I )  /  (  SIGMAX*SIGMAY*PI ) 

T=TSTART 

DELTAT=T/2. 

ILESS=0 

IGREAT=0 

CONTINUE 

IF! ILESS  .EQ.  1  .AND.  IGREAT  .EQ.  1)  THEN 
I LESS =0 
I GREAT =0 

DELTAT=DELTAT/2. 

END  IF 


IF! F(T)  .GT.  .001)  THEN 
T=T-DELTAT 


IGREATOT 
GO  TO  10 

END  IF 


IF(F(T)  .LT.  -.001)  THEN 
TOT-DELTAT 
ILESS= 1 
GO  TO  10 
END  IF 


CINV=T 


RETURN 

END 


FUNCTION  TO  COMPOTE  THE  EXPONENTIAL  POD  FUNCTION. 


FUNCTION  DENSE( I .TIME ) 


COMMON / FACTS / S IGMAX , S IGMAY , CEP .WIDTH ( 2 ) , SPEED i 2) 


PI-3.14159 

C=WIDTH( I ) *SPEED< I ) / < SIGMAX*S IGMAY* PI ) 
ARG = SORT ( COTIME i 
DENSE-C*EXP ( - ARG ) / 2 . 

RETURN 

END 


FUNCTION  TO  COMPUTE  THE  INVERSE  OF  THE  EXPONENTIAL 
POD  FUNCTION. 

X-DENSITY  VALUE  H.E.  THE  ARGUMENT  OF  THE  INVERSE  FUNCTION) 


FUNCTION  DINV(I.X) 


COMMON / FACTS / S IGMAX , S IGMAY . CEP .WIDTH ( 2 ) , SPEED ( 2 ) 


on  on  ooooo 


PI=3. 14159 

C=WIDTH( I ) *SPEED( I ) / ( SIGMAXASIGMAY^PI ) 
ARG=2 . *X/C 

DINV=(LOG(ARG) )A*2/C 

RETURN 

END 


FUNCTION  SUBPROGRAM  TO  SAMPLE  FROM  A  DISCRETE 
DISTRIBUTION. 


FUNCTION  S AMPLE (D 1ST,  ! ARRAY  CONTAINING  FREQ  FUNCTION 

2  J,  (LENGTH  OF  ARRAY 

3  ISEED )  '.SEED  FOR  RANDOM  NUMBER  GEN 


DIMENSION  DIST( J , 2 ) 
COMMON /UN IT / ICRT. ITERM 


Y=RAN( ISEED) 

SUM=0 . 

1  =  1 

100  CONTINUE 

5UM=SUM+DIST  (1,2) 

IF ( SUM  .GE.  Y)  THEN 
SAMPLE-DISK  1,1) 

RETURN 
END  IF 
1  =  1+1 

IF( I  .GT.  J)  THEN 

WRITE( ICRT, A )  'ERROR  IN  SAMPLE' 
ELSE 

GO  TO  100 
END  IF 
END 


$  ty  blnom.for 

OPEN  <  UNIT=4 , NAME= ' IALIVE ' , FORM =' FORMATTED ' ,TYPE='0LD' ) 

READ ( 4 , a ) MMAX 

ISUM=0 

ISUM2=0 

DO  100  1=1, MMAX 
READ( 4 , * ) IPART, IQUIT 
ISUM= ISUM+IPART 


ISUM2=ISUM2+IQUIT 
CONTINUE 
FSUM=FLOAT( ISUM) 

FMAX=FLOAT<MMAX) 

FSUM2=FLOAT( ISUM2 ) 

FRACTION=FSUM/FMAX 

FRACTION 2 =FSUM2/FMAX 

WRITE( 6 ,* ) ' #  OF  SAMPLES  =',MMAX 

WRITE( 6,*) 'FRACTION  FOUND  ALIVE  = ' , FRACTION 

WRITE<6,*) 'FRACTION  PRESUMED  DEAD= ' .FRACTION 2 

ALF=1 . 96 

ROOT=SQRT<  FSUM* ( FMAX-FSUM )  / FMAX*ALF**2 /  4 ) 

TERMl=ALF*ROOT 

TERM2 =FSUM+ALF*  *2/2 

TERM3=FMAX+ALF**2 

BLOWER= ( TERM2-TERM1 ) /TERM3 

BUPPER= (TERM2+TERM1 ) /TERM3 

WRITE<6,*) 'LOWER  BOUND  OF  95%  CONFIDENCE  INTERVAL 
WRITE (6,*) 'UPPER  BOUND  OF  95%  CONFIDENCE  INTERVAL  = 
END 


'  .BLOWER 
'  .BUPPER 


$  lo 


ROEHRIG  logged  out  at  23-APR-1985  15:37 

Remote  disconnect 


